Autoionization of an ultracold Rydberg gas through resonant 
dipole coupling 
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Abstract. We investigate a possible mechanism for the autoionization of ultracold Rydberg gases, based on 
the resonant coupling of Rydberg pair states to the ionization continuum. Unlike an atomic collision where 
the wave functions begin to overlap, the mechanism considered here involves only the long-range dipole 
interaction and is in principle possible in a static system. It is related to the process of intermolecular 
Coulombic decay (ICD). In addition, we include the interaction- induced motion of the atoms and the 
effect of multi-particle systems in this work. We find that the probability for this ionization mechanism 
can be increased in many-particle systems featuring attractive or repulsive van der Waals interactions. 
However, the rates for ionization through resonant dipole coupling are very low. It is thus unlikely that 
this process contributes to the autoionization of Rydberg gases in the form presented here, but it may 
still act as a trigger for secondary ionization processes. As our picture involves only binary interactions, 
it remains to be investigated if collective effects of an ensemble of atoms can significantly influence the 
ionization probability. Nevertheless our calculations may serve as a starting point for the investigation of 
more complex systems, such as the coupling of many pair states proposed in [T]. 

PACS. 32.80.Zb Autoionization - 32.80.Ee Rydberg states - 34.20.Cf Interatomic potentials and forces 
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1 Introduction 

Ultracold Rydberg gases have been found to autoionize, 
typically on time scales of fis [2]. Different processes have 
been identified as possible mechanisms leading to this ion- 
ization. Direct photoionization by black-body radiation is 
one of the important processes. However, under the typical 
conditions of ultracold experiments, the photoionization 
rates are still dominated by interaction-induced collisions 
of Rydberg atoms, which lead to Penning ionization, 



Ryd(n) + Ryd(n) — > Ryd(n') + ion + e 



(1) 



In this process the long-range interactions accelerate the 
initially ultracold Rydberg atoms towards each other until 
they collide, and the ionization cross sections are typically 
very large |3|4j . Recent experiments have identified colli- 
sional ionization to be one of the main ionization mech- 
anisms of ultracold Rydberg gases |5l6j . although other 
processes also seem to be involved in the dynamics of 
Rydberg autoionization. Experiments with dense samples 
show that ions can be observed even before collisions may 
have occurred [Tj. This is attributed to a series of reso- 
nant dipole-coupled pair states which may allow one of 
the atoms of a pair to ionize within a time much shorter 
than the timescale of interaction-induced motion. In more 
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Fig. 1. Resonant coupling to the continuum. Two atoms are 
prepared in state nl at distance R. Due to dipole coupling, one 
atom can be ionized, while the other one is transferred to n'l' 
with n*' < n* /\/2. 



dilute gases, simulations of atom trajectories in many- 
particle systems are found to underestimate the observed 
ionization rate at high principal quantum numbers [7]. 
Apparently there exists a distance-dependent coupling of 
two-particle or many-particle states to other states involv- 
ing the continuum. This dipole coupling induced autoion- 
ization may also be relevant for the stability of Rydberg 
macromolecules |S]. 

Here, we consider the simple case of a direct coupling 
of a pair state to the ionization continuum (see Fig. [I]). 
We restrict ourselves to the range of distances where the 
wavefunctions do not yet overlap. Once the atoms reach 
the overlap region, the ionization probability can be en- 
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hanced by orders of magnitude 9 . Two atoms, both in 
state \nl) and separated in space by a distance R, are 
coupled by the dipole-dipole interaction in such a way 
that one atom is ionized, while the other one undergoes a 
transition to a lower Rydberg state \n'V). Such a process 
becomes possible whenever the energy difference between 
the \nl) and \n'V) states is larger than (or at least equal 
to) the binding energy of the Rydberg electron in the \nl) 
state. The rate of this process depends strongly on the ini- 
tial Rydberg state and on the distance between the atoms. 
The process has been discussed for thermal gases [TO] and 
for dense cold gases [11] . The mechanism is related to the 
interatomic or intermolecular Coulombic decay (ICD), a 
collective autoionization process which has been intensely 
studied in clusters |12|13|14|15j , and has recently been ob- 
served between OH~ and H 2 molecules [TB]. ICD occurs 
when an excited particle relaxes into a lower state, while 
in a radiationless process a neighboring particle picks up 
the energy and is ionized. While in the typical ICD pro- 
cess one of the involved particles is an ion (see references 
above) or atom [17] in an excited state with only one or a 
few empty low-lying energy levels, highly excited Rydberg 
atoms feature a large number of unpopulated states below 
the Rydberg state. In the present work we investigate how 
relevant such a process can be for the autoionization of 
a dilute ultracold Rydberg gas, when interaction-induced 
motion of the atoms and many-particle effects are taken 
into account. 

The paper is organized as follows: In Section [2] the cal- 
culation of the ionization rates based on dipole interaction 
is described. The ionization probability of a pair of atoms 
on an attractive interaction potential is discussed in Sec- 
tion[3] SectionQ] addresses the influence of the surrounding 
atoms in a cloud on this accelerated pair. Finally, autoion- 
ization in a many-body system of repulsively interacting 
atoms is investigated in Section [5] 
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Fig. 2. Rate coefficient F for different initial states nl of ru- 
bidium (in atomic units): ns (dotted line), np (solid line), nd 
(dashed line). The inset shows the same data in a log-log scale 
to visualize the power law scaling. 





10 2 






: I I I I I I I I 




10 1 






10° 




<D 


10" 1 




0B 


10" 2 




o 


10" 3 




0B 
N 


10" 4 




'c 
o 


10" 5 






10" 6 






10" 7 


I I I I I I I 



1 2 3 4 5 6 7 8 9 10 



distance (um) 

Fig. 3. Distance dependence of the rates r (given in s _1 ) for 
different initial states nl: ns (dotted line), np (solid line), nd 
(dashed line) for n = 100 (upper three traces) and n = 60 
(lower three traces). 



2 Calculation of rates 

Consider a pair of atoms, both in a state \nl), at an inter- 
nuclear separation R. The rate for the transition \ nl, nl) — ► 
\n'V , E, I") can be expressed in terms of the dipole matrix 
element of the bound-bound transition and the photoion- 
ization cross section for the bound- free transition [15] . 
Atomic units are used in the following. Averaging over 
a uniform angular distribution, the rate is given by 

r n l n'l'(R) — — — J^\D n l n >l> f , (2) 

mj n i tn 'i'H 

where D n i „/;/ is the reduced matrix element for the elec- 
tric dipole transition between the two bound states \nl) 
and \n'l'), 0J n i, n >i' > denotes the energy difference be- 
tween these states, and api is the photoionization cross 
section for a transition from \nl) to the continuum with 
the same energy difference. 

The total ionization rate must take into account tran- 
sitions to all lower states \n'V) of one atom, which yield 
sufficient energy for the ionization of the other atom. The 



ionization energy of an atom in the state \nl) with the 
effective quantum number n* is equal to its binding en- 
ergy, Ei on = \E n * | = l/(2n* 2 ). This means that ionization 
is possible for transitions to levels n'* with binding ener- 
gies \E n /*\ — l/(2(n'*) 2 ) larger than the threshold value 
\Eth\ — 2-Eion — l/"* 2 - From here we immediately obtain 
the condition for ionizing levels asri'*< n*/V2. The total 
rate can then be obtained by summing over all states n'* 
that satisfy the above ionization condition: 

r nl (R)= J2 r m,wi'(R), (3) 

n'*<n*/V2 
l'=l±l 

The photoionization cross sections are evaluated using 
the semiclassical formula given in [19120121] . 

apiM = i ^n* s ^ ±1 T^ Dll(UJnl ' n ' l ' h W 

where the sum corresponds to the transitions I — ► I" = 
I + Al into the continuum with L c = max{Z, I"}. The co- 
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efficients 
Dai = 



- — <V (x) + AlJ - $ Ail (x) 

£i X V X 



Table 1. Fit parameters for a power law scaling of the rate 
coefficient of the form t — an b . 



(5) 

are expressed via the conventional Airy function and its 
derivative: 

oo 

$ Afi (x) = -1= J d£ cos(x£ + £ 3 /3 + irAn) (6) 



3X 2/3 



AL = Al 



1 

5L r 



(7) 



(8) 



J nl,n'V 



Here Z\^i denotes the difference between the quantum de- 
fects /ii for series I and I" involved into the photoioniza- 
tion: A/i = [jli+ai — fJ-i- For small and large values of the 
argument x the functions <S>Ap, (%) and & AlL (%) can be ap- 
proximated by simple asymptotical expressions [19j . 

The radial matrix elements for the transitions to bound 
states are calculated numerically by integrating the Schro- 
dinger equation for the corresponding energies using the 
Numerov algorithm [32] . For the low- lying levels the wave- 
functions are determined using a model potential fitted to 
one-electron energies [23] . 

In order to eliminate the R dependence of the rate, we 
define a coefficient r n i, such that 



(9) 



The calculated coefficient r n i for rubidium is plotted in 
Fig. [5] as a function of the quantum numbers n and I. The 
oscillatory behavior of the rate coefficient has the follow- 
ing simple explanation. As discussed above, the ioniza- 
tion can occur if the binding energy E n i* of the final level 
nf* exceeds the threshold value E th . The excess energy 
e = \E n /*\ - \E t h\ = l/2n'* 2 - 1/n* 2 is carried away by 
the emitted electron. Since the photoionization cross sec- 
tion decreases rapidly with increasing energy of the pho- 
toelectron e, the transitions to the \n'V) level closest to 
the threshold will give the main contribution to the ion- 
ization rate. The energy spacing between the bound levels 
increase with decreasing n. Therefore for different \nl) lev- 
els the corresponding highest possible ionizing \n'V) levels 
appear at different separations from the threshold Eth, 
leading to different excess energies e. Hence, the ioniza- 
tion rates for different \nl) levels will be larger or smaller, 
depending on how close to the threshold is the first ioniz- 
ing \n'V) level. 

Neglecting the oscillations, the average rate coefficient 
exhibits a power law scaling with the principal quantum 
number, r = an b , as can also be seen from the logarithmic 
representation in Fig. [2] (inset). The best fit parameters a 
and b corresponding to the different quantum numbers I 
are listed in Table Q] Note that the scaling exponent b can 
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0(8) 


0.014 


5.74 


1 (P) 


0.092 


5.52 


2(d) 


0.240 


5.46 



be calculated analytically for the case of hydrogen [24] , 
yielding b — 5.333 independent of /, which is close to the 
fitted value for rubidium with I = 2. 

Figure [3] shows the calculated rates r n i as a function 
of the interatomic distance R for different quantum num- 
bers. For typical interatomic distances and high principal 
quantum numbers the rates are in the mHz range, so that 
the contribution of this ionization process seems negligible 
for experiments with Rydberg gases. Rates above 100 Hz 
are only possible for very high principal quantum num- 
bers (n > 100) at short distances (~ 1/Ltm). Due to the 
excitation blockade observed for high Rydberg states, the 
minimum distance between Rydberg atoms shortly after 
the excitation is limited, typically to a few /im |25l26j . 
Furthermore, the interaction-induced acceleration will not 
allow the atoms to remain at rest on the timescale of an au- 
toionization process. This phenomenon will be discussed 
in the next section. 

Despite the low rates, the process described here can 
still act as a trigger for additional ionization processes. 
Collisions of other Rydberg atoms with the electron and 
ion produced can lead to secondary ionization and ioniza- 
tion avalanches. In some configurations, these avalanches 
are even found to create ultracold plasmas [2"T] . 



3 Ionization probability during motion 

Two interacting Rydberg atoms will experience a force 
given by the gradient of the interaction potential, which 
can be either attractive or repulsive. Consider an atom 
pair on an attractive van der Waals potential with the 
coefficient Cq such that 

r,i do) 



V v dW = 



i? 6 ' 



The atoms will attract each other and collide within sev- 
eral fis for the typical inital distance of a few /xm. Due 
to the motion of the atoms the resonant ionization rates 
become time-dependent, increasing dramatically as atoms 
approach each other. 

In order to investigate the importance of this mecha- 
nism in such a dynamic system, we calculate the ionization 
probability during the motion. Let p(t) be the probability 
for having been ionized at a time t, and V = p(t co u) be 
the probability for ionization before a collision occurs. V 
can be evaluated by integration, 



dp = (l-p(t))r(R(t))dt 

V = 1 - exp (- [ r(R(t)) dt 



(11) 
(12) 
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As long as the exponent in Eq. 
unity, this can be simplified as 



12|) is small compared to 



r(R(t))dt , 



(13) 



where t co u is the time at which the pair distance is so small 
that the electronic wave functions begin to overlap. This 
distance is set to R co ii — 'in* 2 , below which collisional 
ionization occurs [3]. R(t) is evaluated using the classical 
equation of motion on the potential given by Eq. (fT0|) . The 
coefficient C§ is taken from the calculations in Ref. [28] , 

The collision time t co u for atoms with mass M initially 
at a distance Ro may be evaluated analytically assuming 
that the initial velocity of the atoms is zero. In this case 
the collision occurs with zero impact parameter and t co u 
can be estimated by writing down the equation of energy 
conservation with the reduced mass fi = M/2, 



V{R ) = V{R) 



1 



dR 
~dt 



(14) 



and integrating from R 

Rcoll, 



dt 



Ro to the collision radius R 



dR 



v/2/i- 1 (V(R ) - V(R)) 



(15) 



The collision time does not depend strongly on the exact 
value of Rcoii and can be estimated by setting R co u = 0. 
It is also independent of the exact shape of the interaction 
potential at short distances, as the atoms spend most of 
their time at larger separations. The integral can thus be 
solved for a simple van der Waals potential to yield the 
collision time 

fW 



tc 



0.22 



Rn 



(16) 



For large initial distances (i.e. Rq/R co ii ^> 1, V(Rq) = 0), 
the ionization probability can be written analytically as 



V 



4C 6 



Ro 



dRR 



-3 



Reo 



4C 6 2Rl 



(17) 



Fig. 01a) shows the ionization probability V, calcu- 
lated by numerical integration of Eq. (|13[) . as a function 
of the initial distance for pairs of atoms in the states 60D 
(C 6 = lxlO 21 ), 80D (red, C* 6 = 2.4xl0 22 ), and 100D (red, 
Cg = 2.8 x 10 23 ). The ionization probabilities increase 
slightly with increasing initial distance, approaching the 
value found from Eq. [TTJ In order to reach these proba- 
bilities, the observation time must be at least as long as 
the collision time plotted in Fig.[4^b). For a limited obser- 
vation time the rates will decrease again with increasing 
i?o- Despite the higher value of r for increasing principal 
quantum number n, the probability V decreases, as the 
atoms spend only a short time in close proximity to each 
other before they collide. If we let 



v = r(R(o))t coU 



22 \ F 
V C'e Rq 



(18) 
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Fig. 4. Various parameters of ionizing Rydberg pairs GOD- 
GOD (solid lines), 80D-80D (dashed lines), 100D-100D (dotted 
lines) as a function of the initial interatomic separation Ro: (a) 
resonant dipole ionization probability V before the collision; 
(b) collision time t co u; (c) relative increase of the ionization 
probability V/Vo for a moving pair compared to a static pair. 



be the ionization probability for a static pair interacting 
during the same time t co u, then we obtain the relative 
probability V/Vo (see Fig. EJc)). For large Rq we obtain 
from Eqs. ([T7J) and (HHJ) the relation 



V/V ~l.UR 2 JR 2 oU 



(19) 



For typical initial distances Rq the attractive motion in- 
creases the ionization probability by one to two orders of 
magnitude compared to the static pair. 



4 Influence of surrounding atoms 

So far we have considered the simple case of two atoms. In 
a large ensemble of atoms, e.g., an ultracold atom cloud, 
all the surrounding atoms will contribute to the ioniza- 
tion rate. But as the arrangement of atoms is random and 
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not a crystalline structure, there will always be exactly 
one nearest neighbor. Consider an atom with its nearest 
neighbor at a distance i?n. According to Eq. ([9]), the ion- 
ization rate for this pair is given by 



r ■ =L 

J- pair — j^q 



(20) 



To estimate the importance of the surrounding atoms, we 
assume a constant density p and calculate an additional 
rate r s induced by all other atoms except the nearest 
neighbor by integrating from the radius i?n to infinity (as 
the nearest neighbor is at i?n by definition, all other atoms 
must be at larger distances): 



r, = 



Ro 



4?ri? 2 p r{R)dR 



(21) 



Here r(R) = r/R 6 is the ionization rate of a pair at 
distance R. The integral yields 



4tt ~ 1 



(22) 



For a homogeneous density p the average nearest neighbor 
distance can be estimated as a — (5/9V- 1 / 3 [29]. Ass rim- 
ing Rq — a, we can express p in terms of i?n and obtain 



r -*L r-Xr- 

s 3 \9 R* 







0.718 r, 



pair • 



(23) 
(24) 



The ionization rate induced by all the atoms with R > 
a (only pairwise interactions are considered) is therefore 
smaller than the rate induced only by the nearest neigh- 
bor at Rq = a. The ratio r s /r pa i r will be even smaller for 
very close pairs with Ro < a. If the atoms are accelerated 
due to an attractive van der Waals interaction, the 1/i? 6 
scaling of the potential accounts for a much larger acceler- 
ation towards the nearest neighbor compared to any of the 
surrounding atoms. The nearest neighbors will therefore 
approach each other quickly and the influence of the sur- 
rounding atoms will be even smaller compared to a static 
system. 

In order to take the complete distribution of nearest 
neighbors into account, a Monte Carlo simulation has been 
performed, where a large number of atoms are placed ran- 
domly in space. Fig. [5] shows the results for different sear- 
ings r oc 1/R a of an accumlative quantity r for a ran- 
domly picked atom, taking different numbers N of its neig- 
bors into account. The results shown are the total rates of 
100 randomly picked atoms, scaled to the respective value 
of a single neighbor for easier comparison. For a scaling of 
a = 6 clearly only the first neighbor has an effect, while 
for a < 3 the accumulated rate increases quickly. 

Generally, for quantities r scaling as 1/i? 3 the influ- 
ence of the surrounding atoms becomes significant. This 
can be seen from the fact that the integral in Eq. (|21j) 




3 4 5 6 7 
number of partners N 

Fig. 5. Monte-Carlo simulation for the accumulated rates 
S =i V-^» when the N nearest neighbors are taken into ac- 
count and different scalings a = 1, 2, 3, 6. While for a — 6 the 
rate approaches a limit as N — > oo and is mainly determined 
by the first neighbor, the scalings a < 3 lead to a diverging 
rate where all surrounding neighbors have to be considered. 
All values are normalized to the corresponding N — 1 rate. 



then diverges and can only be evaluated with an upper 
limit R m ax ■ Expressing the limits of the integral in terms 
of the average nearest neighbor distance a with Ro = a 
and Rmax = 0a (J3 > 1) the integral yields 



T 3 , s = 2.155 In (3 r< 



3, pair 



(25) 



where Z3 jS and i^pmr now represent arbitrary accumula- 
tive quantities scaling as 1/i? 3 . For (3 — 100 (a realistic 
value considering a magneto-optically trapped cloud), the 
contribution of the surrounding atoms outweighs the near- 
est neighbor by an order of magnitude. Depending on the 
scaling of the quantity r with the interatomic distance, 
the influence of other atoms surrounding a nearest neigh- 
bor pair may or may not be significant. This should be 
kept in mind when estimating energy shifts of an atom in- 
teracting via van der Waals (1/i? 6 ) or dipole-dipole (1/i? 3 ) 
interactions. In the van der Waals case, it may be suf- 
ficient to consider only the nearest neighbor, while the 
dipole-dipole interaction requires taking a larger number 
of atoms into account. Note that in the case of dipole- 
dipole interactions, the angular dependence also has to be 
considered. 

The estimates presented in this section justify the cho- 
sen approach to consider an attractively interacting sam- 
ple as a collection of pairs, and to derive ionization prob- 
abilities in a pair picture as described above, even if the 
atoms are embedded in a larger ensemble. 



5 Repulsive many-particle dynamics 

The derivation of the resonant ionization rates in terms of 
a two-atom picture is only justified if the atoms are pulled 
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Fig. 6. Ionization probability for pair of repelling atoms in 
state 60S (Ce, = — 10 21 au) for different initial distances: (a) 
2/zm, (b) (c) 6/im. The dashed lines show the calculated 

probabilities for two atoms in free space, the solid lines show 
the results of multi-particle calculations for the same pair em- 
bedded in a cloud of interacting atoms. 



together by an attractive force, thereby reducing the influ- 
ence of the surrounding cloud. If the atoms are prepared 
in states exhibiting repulsive interaction, this assumption 
cannot be made. In this case the dynamics of the whole 
sample of atoms must be considered. An isolated pair of 
atoms will be drawn apart immediately, so that the ioniza- 
tion probability is small. In a many-particle system where 
all atoms are repelling each other, each particle may come 
into the vicinity of other atoms several times while mov- 
ing across the cloud [7J. This increases the average time 
spent at close distance and thus increases the ionization 
probability. In this simplified picture, the atoms do not 
collide at all. We thus need to specify the probability for 
resonant dipole ionization as a function of time, V = V(t). 

Figure [S] shows the ionization probability for a pair of 
atoms in the 60S state (with Ce = — 10 21 au). Two isolated 
atoms in free space will be accelerated away from each 
other, and the ionization probability reaches a steady- 
state value. The same nearest-neighbor pair of repelling 



atoms placed in the middle of a large atom cloud has a 
higher ionization probability, as both atoms come close to 
other particles while moving away from each other. The 
probability to ionize one of these two specific atoms in- 
creases steadily in time. In the calculation presented in 
Fig. [6l a density of Rydberg atoms of 1.5 x 10 9 cm~ 3 is 
assumed, and the atoms surrounding the initial pair are 
placed randomly in space. The dynamics of the system 
is then calculated by solving the equations of motion for 
all atoms according to the van der Waals interaction po- 
tentials in the same way as described in Ref. [7J, and the 
results are averaged over 50 runs of the simulation. For 
typical nearest-neighbor distances of about 4 /xm, the sur- 
rounding atoms in this example enhance the probability 
to ionize within 30 fis by roughly a factor of 1.4. How- 
ever, the total ionization probability within this time is 
of the order of 10~ 10 , which is far from being observable 
exp er iment ally. 

6 Conclusion 

We have presented calculations of ionization rates and es- 
timates of many-particle effects for an autoionization pro- 
cess of ultracold Rydberg atoms based on resonant dipole- 
dipole interaction. Our approach involves only binary in- 
teractions and a direct coupling to a continuum state. 
The estimated rates for this process are not sufficient to 
explain the observed autoionization in dense or strongly 
interacting Rydberg gases not caused by collisions |l|7j . 
Taking into account the motion and the influence of a 
many-particle environment, the expected ionization prob- 
ability for this process increases, but it is still too low 
to be observed experimentally. This suggests that other 
interaction- induced ionization processes exist, which con- 
tribute considerably during time scales shorter than the 
collision time of the atoms. Furthermore, even a process 
with a low intrinsic rate may act as a trigger for secondary 
ionization. 

We believe that collective effects of larger ensembles 
of interacting atoms may have a great influence on dipole 
coupling induced autoionization. As a next step, the model 
could be extended to a coupling of many intermediate pair 
states as proposed by pQ , and for further analysis the sys- 
tem could be described in terms of multi-particle states 
instead of pair states. 

The discussion of the dynamics in many-body systems 
presented here is not restricted to the autoionization pro- 
cess but may also be helpful for other dynamical processes 
involving Rydberg interactions. 
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